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Abstract 

With the help of the Drinfeld twist or factorizing F-matrix for the eight-vertex 
SOS model, we obtain the explicit determinant expression of the partition function 
of the eight-vertex model with a generic non-diagonal reflecting end and domain wall 
boundary condition. Our result shows that, contrary to the eight-vertex model without 
a reflection end, the partition function can be expressed as a single determinant. 
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1 Introduction 



The domain wall (DW) boundary condition of a statistical model on a finite two-dimensional 
lattice was first introduced in |T] for the six-vertex model. The partition function of the 
model (or DW partition function) was then expressed in terms of a determinant [21 [3]. 
Such a determinant representation of the partition function has played an important role 
in constructing norms of Bethe states, correction functions [H |5l [6] and thermodynamical 
properties of the six- vertex model [TJ, and also in analyzing the Toda theories |8J. Moreover, 
it has been proven to be very useful in solving some pure mathematical problems, such as 
the problem of alternating sign matrices [9]. Recently, the determinant representations of 
the DW partition function have been obtained for various models [ini IHl [El [ISl Ull UHl HE] . 

For a two-dimensional statistical model with a reflection end [T7], in addition to the 
local interaction vertex, a reflecting matrix or K-matrix [TS] which describes the boundary 
interactions needs to be introduced at the reflection end of the lattice (see figure 4 below). 
The DW partition function of the six- vertex model with a diagonal reflection end was exactly 
calculated and expressed in terms some determinant Such determinant representation 
of the model was re-derived [20] by using the Drinfeld twist or factorizing F-matrix [2T] 
method f22]. However, it is highly non-trivial to generalize this result to models with non- 
diagonal boundary terms [23l [23 [251 [261 [271 [2HI [23 [301 [311 [321 [33 

Very recently, determinant representations of DW partition functions have been obtained 
for the trigonometric SOS model with reflection end [H] and the six-vertex model with a 
non-diagonal reflection end [12]. 

Among solvable models, elliptic ones stand out as a particularly important class due to the 
fact that most trigonometric and rational models can be obtained from them by certain limits. 
In this paper, we focus on the most fundamental elliptic model — the eight-vertex model 
[l3] whose trigonometric limit gives the six-vertex model. Here we obtain the determinant 
representation of the DW partition function of the eight-vertex model with a non-diagonal 
reflection end specified by the generic non-diagonal K-matrix [^ [15] . Our result shows 
that contrary to the eight-vertex model without reflection end whose DW partition function 
can only be expressed as a sum of determinants [131 [IH [IS] > the DW partition function of 
the eight-vertex model with a non-diagonal reflection end can be represented as a single 
determinant. Such a single determinant expression will be essential for the study of the 
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Bethe vectors of the open XYZ chain with non-diagonal boundary terms. 

The paper is organized as follows. In section 2, after introducing our notation and some 
basic ingredients, we construct the four boundary states which specify the DW boundary 
condition of the eight- vertex model with a non-diagonal reflection end. In section 3, using the 
vertex-face correspondence relation we express the DW boundary partition function in terms 
of the matrix element of the product of the (face type) pseudo-particle creation operators. 
In section 4, with help of the F-matrix of the eight-vertex SOS model we give the completely 
symmetric and polarization free representations of the pseudo-particle creation operators in 
the F-basis. In section 5, we obtain the determinant representation of the DW partition 
function. In section 6, we summarize our results and give some discussions. 

2 Eight-vertex model with a reflecting end 

In this section, we briefly review the DW boundary condition for the eight- vertex model with 
non-diagonal reflecting end on an iV x 2N rectangular lattice. 



2.1 The eight-vertex R- matrix and associated K- matrix 

Let us fix T such that Im(r) > and A,- be the lattice generated by 1 and r. Introduce the 
following elliptic functions 
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(m,2t), j = 1,2; 

The cr-functioiij satisfies the so-called Riemann identity: 

a{u + x)(y{u — x)a{v + y)a{v — y) — (y{u + y)a{u — y)a{v + x)a{v — x) 
= a{u + v)a{u — v)a{x + y)cr{x — y), 



(2.1) 
(2.2) 



(2.3) 



which will be useful in the following. Moreover, for any a = (0:1,02), 01,02 ^ ^2, we can 
introduce a function a„(n) as follow 
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"'^Our (7- function is the z?-function i9i(w) It has the fohowing relation with the Weierstrassian a- 

function au]{u): cr^iu) cx e^^^^a{u) with 771 = •7r^(^ — 4^^-,^ T^q^) ^"^^ 1 — ■ 



The above definition implies the identification (T(o,o)(w) = <y{u). 

Let \^ be a two-dimensional vector space and {ei\i = 1, 2} be the orthonormal basis of 
V such that (ej, ej) = 6ij. The well-known eight- vertex model R-matrix R{u) G End(V^(8> V) 
is given by 



R{u) 



( a{u) d{u) \ 

h{u) c{u) 
c{u) h{u) 
\ d{u) a{u) j 



(2.5) 



The non-vanishing matrix elements are 



_ 0(1) {u) 0(°) {u + T]) a{T]) _ 0(0) {u) 0(1) (u + J]) a[r]) 

''^''^ ~ 0(i)(O)0(o)(r/)a(n + r7)' ~ 0{i)(O) 0(o)(r/) a(u + 77) ' 

^ 0(^)(t^)0(^)(tx + 77)^(77) ^ g(o)(^)g(o)(^ + ^)^(^) 

0(i)(O)0(i)(r/)a(M + ?7)' 0(i)(O)0(i)(r/)(T(« + r7)- 

Here u is the spectral parameter and j] is the so-called crossing parameter. The R-matrix 
satisfies the quantum Yang-Baxter equation (QYBE) 

^1,2(^1 - U'^Rx^ziux - Us)R2,3iu2 - U3) = i?2,3(M2 " U3)Ri^s{Ui - U3)Ri^2{Ul " U2) . (2.7) 

Throughout we adopt the standard notation: for any matrix A G End(y), Aj (or yl-')is an 
embedding operator in the tensor space V ® ■ ■ ■, which acts as A on the j-th space and 
as identity on the other factor spaces; Rij (u) is an embedding operator of R-matrix in the 
tensor space, which acts as identity on the factor spaces except for the i-th and j-th ones. 

For a model with reflection end |18], in addition to the R-matrix, one needs to introduce 
K-matrix K{u) which satisfies the reflection equation |17] (RE) 

Rl,2{ui - U2)Ki{ui)R2,l{Ui + U2)K2{U2) 

= K2{U2)R1,2{U1 + U2)Ki{Ui)R2,l{Ui - U2). (2.8) 

In this paper, we consider the K-matrix K{u) which is a generic solution H5] to the 
RE (12. 8p associated with the R-matrix (12. 5p . 

K{u) = ko{u) + K{u) a^' + kyay + h{u) a^ (2.9) 

where a^,ay,a^ are the Pauli matrices and the coefficient functions are 

ko{u) 



a{2u) a(Ai + A2 - |) a(Ai + C) ^(As + C) 



2 a{u) a{-u + Ai + A2 - |) a(Ai + C + u) ^(^2 + C + «) ' 
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kxiu) 

ky{u) 

kz{u) 



a{2u) o-(i,o)(Ai + A2 



(1,0) (Ai + C) ^(1,0) (A2 + C) 



2 a(i,o)(n) a{-u + Ai + A2 - |) a(Ai + C + u) (^{>^2 + C + u)' 

i a{2u) a(i,i)(Ai + A2 - |) cr(i^i)(Ai + C) q'(i,i)(A2 + C) 
2 cr(i,i)(M) a{-u + Ai + A2 - |) a{Xi + ( + u) a{\2 + C + u)'' 

(t{2u) a(o,i)(Ai + A2 - I) o-(o,i)(Ai + C) cr(o,i)(A2 + C) 



(2.10) 



2 cT(o,i)(n) cr{-u + Ai + A2 - I) a(Ai + C + ^(A2 + C + u)' 
Here there are 3 free boundary parameters {Ai, A2, C}) which are related to different boundary 
interactions. It is very convenient to introduce a vector X E V associated with the boundary 
parameters {Aj}, 

2 

X = Y,^kek. (2.11) 



k=l 



2.2 The model 

The partition function of a statistical model on a two-dimensional lattice is defined by 



Z 



E 



where E is the energy of the system, k is the Boltzmann constant, T is the temperature 
of the system, and the summation is taken over all possible configurations under certain 
boundary condition such as the DW boundary condition. The model we consider here has 
eight allowed bulk vertex configurations 
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Figure 1. Vertex configurations and tlieir associated Boltzmann weights, 
and four allowed configurations at each reflection end 
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Figure 2. Reflection ends and the associated Boltzmann weights. 



where 1 and 2 respectively denote the spin up and down states. Each of the eight bulk 
configurations is assigned a statistical weight (or Boltzmann weight) Wi, while each of the 
four reflection configurations is assigned a weight bi. Then the partition function of the 
model with a reflection end can be rewritten as 

Z = J2 ^'i"' ^2"' w's"' w^4"' w^s"' We""" wr""' u^s"' b[' b^^ b^i 

where the summation is over all possible configurations with and Ij being the number 
of vertices of type i and the number of reflection ends of type j respectively. The bulk 
Boltzmann weights which we consider here have Z2-symmetry, i.e., 

a = wi = W2, b = W3 = W4, c = W5 = wq, d = wi = ws, (2-12) 

and variables a, b, c, d satisfy a function relation, or equivalently, the local Boltzmann weights 
{wi\ can be parameterized by the matrix elements of the eight-vertex R-matrix (12. 5p as in 
figure 3. At the same time, the weights {bi} corresponding to the reflection end can be 
parameterized by the matrix elements of the corresponding K-matrix (12. 9p as in figure 3. 
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Figure 3. The Boltzmann weights and elements of the eight-vertex R-matrix and K-matrix. 

Then the corresponding model is called the eight-vertex model with a reflection end. Its 
partition function is thus given by 



As shown in figure 3, to parameterize the bulk Boltzmann weights in terms of the elements 
of the R-matrix, one needs to assign spectral parameters u and ^ respectively to the vertical 
and horizontal lines of each vertex of the lattice. In an inhomogeneous model, the statistical 
weights are site-dependent. Hence two sets of spectral parameters {uoi} and {fj} are needed, 
see figure 4. The horizontal lines are enumerated by indices 1, . . . , with spectral parameters 
{fj}, while the vertical lines are enumerated by indices 1, . . . , with spectral parameters 
{ua} (The 2N parameters {"Uq} are assigned as follow: U2i = Ui and "^24+1 = —Ui, as shown 
in figure 4.). The DW boundary condition is specified by four boundary states |f2'^^)(A)), 
|0(i)(A)), (fi(^)(A)| and (n(2)(A)| (the definitions of the boundary states will be given later, 
see fl2.38p - (l2.4ip below). These four states correspond to the particular choices of spin states 
on the four boundaries of the lattice . 
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Figure 4. The eight-vertex model with a non-diagonal reflection end and the DW condition. 

Some remarks are in order. The boundary states not only depend on the spectral param- 
eters (I^^^^^A)) and (^^(^)(A)| depend on {^i}, while and (0(2)(A)| depend on 
but also on two continuous parameters Ai and A2. However, in the trigonometric limit (i.e., 
setting A2 = I and then taking r — j- +ioo), the corresponding boundary states |r2*^-'^)(A)) and 
(n'^"'^^(A)| (or |n^^^(A)) and (fi*^^^(A)|) become the state of all spin up and its dual (or the state 
of all spin down and its dual) up to some over-all scalar factors. 

The partition function of the eight- vertex model with a non-diagonal reflection end spec- 
ified by the generic K-matrix K{u) ( \2.9\\ and the DW boundary condition is a function of 
2N + 3 variables {mq,}, {^j, Ai, A2 and (, which is denoted by Z^{{ua}', A; Q. Due to 
the fact that the local Boltzmann weights of each vertex and reflection end of the lattice are 
given by the matrix elements of the eight-vertex R-matrix and the associated K-matrix (see 
figure 3), the partition function can be expressed in terms of the product of the R-matrices, 
the K-matrices and the four boundary states 



^iv(K};{6};A;C) = 
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xRi^n{ui-^n) ■ ■ ■ Ri,i{ui-^i)Ki{ui)Rij{ui+(i) . . . Rn,i{ui+^n) 

Xi?^ jv(m7V-C7v) • • • RN,l{uN-^l)Kf^{uN)Rl,N{uN+^l) ■ ■ ■ Rn,n{'^N+^n) 

x|n(i)(A))|fi(2)(A)). (2.13) 

One can rearrange the product of the R-matrices in (I2.13P in terms of a product of the 
so-called double-row monodromy matrices 

Zn{M, m A; 0={^^'\m^^'\^)\ Mui) . . . |0«(A))|r](2)(A)), (2.14) 

where the monodromy matrix Ti{u) is given by 

Ti{u) = Ti{u;^i, . . . ,^n;C) = RiN{ui-^N) ■ ■ ■ Ri,i{ui-^i) 

xKi{ui)Ri-i{ui+Ci) . . .RN,i{ui+CN)- (2.15) 

The double-row matrix T{u) has played an important role in constructing the transfer matrix 
for an open spin chain [18]. The QYBE fl2.7p of the R-matrix and the reflection equation 
(12. 8p of the K-matrix ensure that the monodromy matrix Ti{u) satisfy the following exchange 
relation 

Rijiui - Uj) Tj{ui) Rj-i{ui + Uj) Tj{uj) = Tj{uj) Rijiui + Uj) Tj{ui) Rj-i{ui - Uj). (2.16) 

2.3 The boundary states 

From the orthonormal basis {ej} of V, we define 

1 ^ ^ 
i = ei-e, e=2X]^'^' ^ = ^'2) then^? = 0. (2.17) 

k=l i=l 

Let f) be the Cartan subalgebra of Ai and f)* be its dual. A finite-dimensional diagonalizable 
f)-module is a complex finite-dimensional vector space W with a weight decomposition W = 
©^g(,. so that I) acts on W[fi] by xv = fi{x)v, {x E i), v E W[fj]). For example, the 
non-zero weight spaces of the fundamental representation Va^ = = V are 

W[i\ = Cei, i = 1,2. (2.18) 

For a generic m E V, define 

nii = (m, ei), niij = rui - nij = (m, - e^), ij = 1, 2. (2.19) 

9 



Let R{u, m) G End(y ® V) be the R-matrix of the eight- vertex SOS model [13] given by 

2 2 

R{u- m) = Y^ R{u- m)lEu®Eii+Y^ {^(«; m^Eii^Ejj+R{u; m)i^E,i^Eij} , (2.20) 

where Eij is the matrix with elements {Eij)[ = SjkSu. The coefficient functions are 

Riu;m)t = l, R{u; m):j = " ^) , , ^ ^, (2.21) 

+ 77) o-(mij) 



+ (2.22) 

o-(m + 7]) a[mij) 



and is defined in fl2.19p . The R-matrix satisfies the dynamical (modified) quantum 
Yang-Baxter equation (or the star-triangle relation) 



Ri^2{ui -U2]m- r]h^^^)Ri^s{ui - %; m)i?2,3(^^2 -Us;m- r]h^^^) 

= R2,z{u2 - U3] m)Ri^s{ui -U3;m - r]h^'^^)Ri^2{ui - U2] m). (2.23) 

Here we have adopted the convention 

i?i,2(M, m - r]h^^^) vi ^ V2 ^ V3 = {R{u, m-r]^)® id) Vi®V2® V3, if V3 G W[fi]. (2.24) 

Moreover, one may check that the R-matrix satisfies the weight conservation condition, 

[h^^^ + h^^\ Ri^2{u;m)] =0, (2.25) 

the unitary condition, 

Ri,2{u; m) R2,i{—u; m) = id (g) id, (2.26) 

and the crossing relation 

R{u, - El e, + ^) ^(-^ m)u > (2-27) 

where 

= 1, £2 = -1, and 1 = 2, 2 = 1. (2.28) 

Let us introduce two intertwiners which are 2-component column vectors <f)m,m-rij{u) 
labelled by 1, 2. The A;-th element of (j)rn,m-r^j{u) is given by 

€]m-,3{u) = e^'\u + 2m,l (2.29) 
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^W(M + 2m2) 
^(2)(^i + 2m2) 



(2.30) 



where the functions 0^^\u) are given in fl2.2p . Exphcitly, 

/ 6'W(M + 2mi) 
\ e^^>{u + 2mi) 

One can prove the following identity 

^W(n + 2mi) ^W(M + 2m2) 
^(2)(^^ + 2mi) 6'(2)(ix + 2m2) 

where C(t) is non- vanishing constant which depends on r. This implies that the two inter- 
twiner vectors (t)m,m-r,i{u) are linearly independent ioi a generic m&V. 

Using the intertwiner vectors, one can derive the following face-vertex correspondence 
relation |13] 



det 



C(r) a{u + mi + m2- -) cr{mi2), 



(2.31) 



Rl^2{Ul - U2)(pln,m-'ni(^l)'Pm-rii,rn-r^ii+3)i^2) 



fc^(^i)C™,-„/:(«2). 



(2.32) 



A;,/ 



Then the QYBE (12. 7p of the vertex-type R-matrix R{u) is equivalent to the dynamical Yang- 
Baxter equation (I2.23P of the SOS R-matrix R{u, m). For a generic m, we can introduce other 
types of intertwiners 0, (f) which are both row vectors and satisfy the following conditions, 



(2.33) 



from which one can derive the relations, 

2 



[u) = id, 



(2.34) 



-u) = id. 



(2.35) 



Following the method in [48], we check that the K-matrices K{u) given by (12. 9 p can be 
expressed in terms of the intertwiners and diagonal matrices ]C{X\u) as follows 



(2.36) 



Here the diagonal matrix }C{X\u) is given by 



/C(A|«) ^ Diag(MA|n)i, fc(A|n)2) = Diag( ^|^^ + ^ , ^T^^^f-^)- (2.37) 
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Although the vertex type K- matrix (u) given by (12. 9 p is generally non-diagonal, after the 
face- vertex transformation (I2.36p . the face type counterpart IC{X\u) becomes diagonal. This 
fact enabled the authors in [261 EH] to diagonalize the transfer matrix of the open chains with 
non-diagonal terms by applying the generalized algebraic Bethe ansatz method developed in 
[28]. 

Now we are in the position to construct the boundary states specifying the DW boundary 
condition of the eight- vertex model with a non-diagonal reflection end, see figure 4. For any 
vector m e y, we introduce four states which live in the two N-tensor spaces of V (one is 
indexed by 1, . . . , and the other is indexed by 1, . . . , iV) or their dual spaces as follows: 

|(](2)(m)) = </>^^„_^2te)C-^2,rn-2,72(^2)...</>^^_^(^_l)2,„_^^2(^^)' (2-38) 
^m~iN~2)rii,m~(N-2)r}i-rii^~'^'^^ ^m-iN-4:)rii,m-{N-4)r}i~rii^^'^'^') 

X • • • C+7Vr,i,m+7v^i-^i(-^^)' (2-39) 



{Q^^\m)\ = 0;„^^_^i(6)0L^i,„_2^i(6)---0^_^(;v-i)i„M(^^)' (2-40) 

X • • • <^m+r,(Af-2)i,m-r,(Af-2)i-,?2(^^)- (2-41) 

The boundary states corresponding to the DW boundary condition can be obtained through 
the above states by special choices of m and i (for example, m is specified to A which is related 
to the parameters of the K-matrix K{u)). The DW partition function ZN{{ua}; {^i}] A; Q 
(I2.13P of the eight-vertex model with a non-diagonal reflection end becomes 

ZN{,{ua}] {6}; A; C) = 

K) . . . 



xR^n{un - Cn) ■■■ Rn,i{un - 6) K^{un) Ri^n{un + 6) • • ■ Rn,n{un + Cn) 

X'^A,A-r72(?l) • • • ^A-(Af-l)r)2,A-Af»72(?^) X-{N-2)iji,X-{N~l)'ni(~'^'^ ' ' ' ^ X+Nrii,X+{N-l)rii(~''^^) ' 

(2.42) 
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3 Partition function in terms of the face type mon- 
odromy matrix 



Let us introduce the face type one-row monodromy matrix 
Tpil\u) ^ To^i...^(/|n) 



N-l 



Ro,n{u -^n;1-vJ2 ^^'^) • • • ^o.2(^^ - 6; ^ - vh^^^)Ro,i{u - 6; 0, 

i=l 

T^ml TFil\u)l\ 

TFil\u)j Tj.il\u)l J ^^-'^ 
where I is a generic vector in V. The monodromy matrix satisfies the face type quadratic 
exchange relation |l9l [50]. Applying rp(/|'u)* to an arbitrary vector l^i, . . . ,1^) in the N- 
tensor product space V^^ given by 

\tu...,iN) = ei...el, (3.2) 

we have 

TF{l\uyAii, . . . ,iN) = TF{'m]l\u))\ii, . . . .In) 



N-l 

ojv_i...ai k=l 

xi?(^.-6;/-r/?X|i?(«-ei;0r^ (3.3) 

where m = I - rj Ylk=i ^k- 

Now we compute the partition function Zjq{{ua\; {ii}] A; Q (I2.13p . The expression (12.421) 
implies that 

Zn{{Uo}] Hi}; A; C)={^^^\^)\ '^L(^_l)^i+r,2,A-(7V-l)r,i(«l) '^L(^-2)^i,A-(^-l)^i ("^O 

(3.4) 

With the help of the crossing relation (I2.27p . the face- vertex correspondence relation (I2.32p 
and the relations (I2.33p . following the procedure in [28| |5T] . we find that the partition func- 
tion Z7v({uq,}; {^i}; A; Q can be expressed in terms of the face-type double-row monodromy 
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operators 

ZN{{ua}\ {6}; A; C) = (1, . . . , l|r/(A - 2(M - l)r/i, \\ui)l . . . Tf{\, \\um)1 

xTf{\ + 2r]l X\uM+i)l . . . TpiX + Nr]i, X\un)1\2, . . . , 2). (3.5) 

The above double-row monodromy matrix operator {m, X\u)'l is given by the one-row 
monodromy matrix operator Ti?(m; [52j 



^(^2l) f-J^ (^{U + Ck + V) 

<^[Xi + C, + u) 



a{X2 + C + u 

F-basis 



-Tp(m + 2r]2, A|M)^Tp(m + 77!, A + r^l] - n - 77)! } . (3.6) 



In this section, after briefly reviewing the result ^3] about the Drinfeld twist [21] (factorizing 
F-matrix) of the eight- vertex SOS model, we obtain the explicit expression of the double rows 
monodromy operator Tf{m, X\u)1 given by (13.61) in the F-basis provided by the F-matrix. 

4.1 Factorizing Drinfeld twist F 

Let Sn be the permutation group over indices 1, . . . , and = 1, . . . , — 1} be the set 
of elementary permutations in S^. For each elementary permutation Sj, we introduce the 
associated operator ^ on the quantum space 

i-l 

Rl ^{l) ^ R^'il) = i?,,+i(C, - - 77^^ /i(^)), (4.1) 

fc=i 

where / is a generic vector in V. For any s, s' G Sn, operator Rf^' associated with ss' 
satisfies the following composition law: 

Let s be decomposed in a minimal way in terms of elementary permutations, 

s = sp,...sp^, (4.3) 
14 



where /3i = 1, . . . , N — 1 and the positive integer p is the length of s. The composition law 
(14. 2p enables one to obtain operator associated with each s G Sjy. The dynamical 

quantum Yang-Baxter equation (12.231) . weight conservation condition fl2.25p and unitary 
condition (I2.26p guarantee the uniqueness of -Rf jv- Moreover, one may check that 
satisfies the following exchange relation with the face type one-row monodromy matrix (13.11) 

RL.Nmo^,i...Nm = T,%,..N)mRi.Ai - vh^'^), vs e s^. (4.4) 



Now, we construct the face-type Drinfeld twist Fi j^{l) = Fi jv(/; ^i, . . . , ^n) □ on the 
iV-fold tensor product space y®^, which satisfies the following three properties: 

I. lower — triangularity; (4-5) 

II. non — degeneracy; (4-6) 

III. factorizing property : RIj^{1) = F;^^Ij^^{1)Fi,„n{1), Vs G 5^. (4.7) 

Substituting (14. 7p into the exchange relation (14. 4p yields the following relation 

F;l^^^{l)F,MlKi...Nm = T,%...^^{l\u)F;l ^^{l - vh^'^)F,,Ml - Vh^'^). (4.8) 
Equivalently, 

F'iMlK,i...NilW)F,-}Al - Vh'^'^) = i^s(i...7V)(0rJ:.(i...^)(%)i".a...;v)(^ - Vh^''>)- (4.9) 
Let us introduce the twisted monodromy matrix Tq^ j^{l\u) by 

fo^,i...Nm = F^MiK,i...NmFr:Ai'vh'^'^) 

Then (14. 9 p implies that the twisted monodromy matrix is symmetric under Sn, namely, 

fo^,i...Nm = %ii...N)m, yseSN. (4.11) 

Define the F-matrix: 



N 



se5jv {aj}=l i=i 

^In this paper, we adopt the convention: Fs(i...Ar)(0 = -^s(i...Af)(^; Cs(i) j • ■ • : £.s{N))- 
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where is the embedding of the project operator Pa in the i*^ space with matrix elements 
{Pa)ki = ^ki^ka- The sum Y^* in (I4.12p is over all non- decreasing sequences of the labels as{iy. 

> oisii) if s{i + 1) > s{i), 

asii+i) > asii) if s{i + 1) < s{i). (4.13) 

From f l4.13p . Fi...Ar(0 obviously is a lower-triangular matrix. Moreover, the F-matrix is non- 
degenerate because all its diagonal elements are non-zero. It was shown in [53j that the 
F-matrix also satisfies the factorizing property (14. 7p . 

4.2 Completely symmetric representations 

In the F-basis provided by the F-matrix ( I4.12p . the twisted operators Tp{l\u)l defined by 
(I4.10p become polarization free [531 El] ■ Here we present the results relevant for our purpose 



o-("-gt) 



TF{l\u)i = —j ^ \ „ , , ^ , (4.14) 

\ / {«) 

f,m\ = Y.-¥^f-^^ ^^-^^-^^^-^^ , ■ (4.15) 



(i) 



Applying the above operators to the arbitrary state |zi, . . . ,ZAr) given by (13. 2 p leads to 

fp{mAu)\= . -{"-«^+'^) ^ 1 , (4.16) 

o {h - mx - f]) \ 1 ' 



^ a(^x - + ?7)a(/i2) 



x^is^j^il . (4.17) 

^ / 0) 

It then follows that the pseudo-particle creation operator (13. 6p in the F-basis has the following 
completely symmetric polarization free form: 

Tp (m, \\u)\ = — — — 

(r{mi-\2)l}^a{u + ik + ri) 



^ a{Xi + C + u)a{X2 + C + u)a{u - + T])a{u + ^j) 
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o-(«-gj )<T{u+^j +r])a{^i -^j+v) 

xEl2^j^^\ ^(^-^j+^M^+^jM^i-^j) j . (4.18) 

^ iJ) 

5 Determinant representation of the partition function 

In this section we compute the DW partition function ZN{{ua}] {Ci}] ^] C) (13. 5p using the 
expansion of the twisted operator Tfijn, \ \u)\ (I4.18p . 

5.1 Symmetric expression of the partition function 

From the definition of the F-matrix -Fi...Ar(/) (I4.12p . we can show that the state |2, . . . , 2) and 
the dual state (1, . . . , 1| are invariant under the action of Fi jv(/), namely, 

Fi...^(/)|2,...,2) = |2,...,2), (5.1) 
(l,...,l|Fi...jv(/) = (1,...,1|. (5.2) 

Hence the DW partition function ZN{{ua}] {^j}] A; Q can be expressed in terms of the twisted 
operator Tf{m, X\u)l as follow 

ZN{{uah {6}; A; =(1, • • • , l\rpiX-iN-2)rji, \\ui)l . . . TfiX+Nr]!, \\un)1\2, . . . , 2) 

= (1, . . . , 1| Fi...^(A - Nrji) Tf{\ -{N- 2)7^1, \\ui)l . . . 
xTp{\ + iVr/i, X\un)l F-l^{X + Nril) |2, . . . , 2) 

= (1, . . . , l\ff{X-{N-2)ril, X\ui)l . . . ffiX+Nrii, X\un)1\2, . . . , 2). 

Substituting the polarization free expression fl4.18p of the twisted operator Tfi^m, X\u)l into 
the above equation, we have 

ZMW.}: {&}; A; c) = fi "a '"■.1:ta""f''/f U U T'^lt!' ' ' ' • 'I 



xf ^ ^(Ai+C-^.).(A.+C+e.).(2n,).(,) Ey, J^^^S^^^:f ^ 
^a(\i+C+n])cT(A2+C+ni)a(7Xi-ei+r?)a(ni+e,) ^^^^^'^ 1^ 



A a(Ai+C-e*)fT(A2+C+e*)(7(2n^)a(r/) . /a(r.^g>(^,^g,+^V(g,^g,+»7) ^ 

x|2,...,2). 
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Expanding the sum in the above equation corresponding to the spectral parameter un yields 



k: 



\ cr(Ai2 + kri)a{Xi2 - kr] + r]) j-^ jj^ a{ui + + v) 



q-(Ai+C-gi)q-(A2+C+^0cr(2nAf)o-(T/) a{ui - ii)cr{ui + + y) 

^a^i+C+UN)(T{X2+C+UN)(T(uN-Ci+'n)<^(^N+^^ fj^ cr{ui - + r))a{ui + ^i) 



x|2,...,2). 

Iterating the above procedure, we obtain the complete symmetric expression of the partition 
function ZN{{ua}; {6}; A; C) 

7 (I. 1 Jf 1 \ r\ T\ ^(^12 + 2fc?7)a(Ai2 - ^kr] + r]) (T{ui + ^i) 
^iv(IM,l6l,A,C) - 11 ^M.„4-i.^wa.._t..^^^ 1111 



IJl fx(Ai2 + kr])a{Xi2 -krj + T]) ^ 1 11 + + 

xZ^(K};te};A;C), (5.3) 
where the normalized partition function Zj^{{ua}; {Ci}', A; C) is 



7 1 r.i ^n-\"TT^ a(Ai+C-ei.(„XA2+C+eM„X2t^n)^(^) 

5.2 Recursive relation and the determinant representation 

From the expression fl5.4p . it is easy to check that the partition function Ziy{{ua}', A; () 
is a symmetric function of {ua} and {^i} separatively. Moreover, we can show that the 
partition function ZN{{ua}', {Ci}] A; C) satisfy the following recursive relation 

7 ({ \ if^l \ n = tT(Ai+C-6)^(A2+C+6)^(2Mjv)ff(r/) 
^ivU^.I, K.I, A, y ^ ^(A,+^+^^)a(A2 + C+«^)cr(«^-6+^)^(«^+e.) 
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^ TT (^{Ui - ii)a{ui + 6 + ^) TT 0^(0 - 6 + ^) 

xZN-i{{ua}a^N; A; ()• (5-5) 

One can show that the initial condition Zo{{ua}', {^i}] -^S C) = 1 ^"^^ the recursive relation 
(I5.5P uniquely determinate the partition function Z^{{ua}', C) for any positive integer 

N. This fact allows us to obtain the following determinant representation of the normalized 
partition function ZN{{ua}', {6}; 0- 

7 (f 1 Jtl \ A^ _ na=i n!Ii ^(^g - 6)^(^a + ^» + ^) detA/"({na}; {^j}) 



where the N x N matrix //{{ua}] {^i}) is given by 



a(r^)a(Ai + C - Q 

a{Ua - Qcr{Uo, + + r?)a(Ai + C + Ua) 



We shall prove the above determinant representation as follows. Firstly let us intro- 
duce two sets of functions {Bi{{ua}', {^i}', A; C) I = 1? • • • > ^} and {Fj{{ua}', {^i}', A; C) I = 
1, . . . , A^} defined respectively by 

Fi{{ua}; {6}; A; C) 



a(r/) 



xdet 



Then the proof of (15. 6p is equivalent to proving the equation 



(5.9) 



{6}; A; C) = Fi{{ua};{Ci};>^;0, for any positive integer I. (5.10) 
We prove (I5.10p by induction. 

• Prom direct calculation, we can show that (I5.10p holds for the case of / = 1, namely, 



5i(Mi;ei;A;C) = FiK;ei;A;C) 
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Suppose that (15.101) holds for the case of / < A^— 1. We now prove that it is satisfied also 
for I = N. It is easy to check that both Bj^{{ua}] {^i}', A; () and Fj^{{ua}] {^i}', A; () 
are symmetric functions of {ua}- Hence it is sufficient to prove that they are equal 
as functions of uj^. The recursive relation (15. 5p of Z]\r{{ua}', {^i}] A; C) implies that 
BnHuo}', {^i}', A; C) satisfies the following relation 

X n ^^17^ B^^.iiu^U,,; A; C). (5.11) 

The determinant representation of the function Fiy{{ua}', A; C) implies that it sat- 
isfies the following recursive relation 

-^(T(MAr-4i + r;)a(MAr + 4i) a{UN-Ui)a{UN+Ui + Tj) 

j4 f^(4i-4i)c^(4j+6) 

The determinant representation (15. 9p and its recursive relation (I5.12p of F^{{ua}', A; C) 
and the recursive relation (15. lip of BM{{ua}] {6}; A; C) imply that these two functions, 
as functions of un, have the same simple poles located at|j: 

6-7], -6 mod(A.), i = l,...,N. (5.13) 

Direct calculation shows that the residues of the two functions at each simple pole 
(I5.13P are indeed the same 

Res„^=^ {FNi{ua}; {6}; A; C))=Res„^=^ {BNi{ua}; {6}; A; C)) , z=^~r],-^i. (5.14) 

Let us introduce the function f{u): 

f{uM) = Fjv(K}; {6}; A; C) - B^iM; {6}; A; 0- (5.15) 

The quasi-periodic properties of the cx-function, 

a{u + l) = -(x{u), (T(M + r) = -e-2^"("+5)a(M), (5.16) 



■^The determinant expression (|5.9p guarantees that the apparent poles in (|5.12p . which are located at 
ui, —ui — rj mod(AT-) for I = 1, . . . , iV — 1, do not really be poles. 
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allow us to derive the following quasi-periodic properties of the function f{u) 

f{u + 1) = f{u), f{u + r) = e-2^-(-2"-^--)/(w). (5.17) 

This indicates ^43j that the number of zero points minus that of poles of the function 
f{u) inside the fundamental domain generated by 1 and r should be —2. However fl5.13p 
and (15.141) imply that the function f{u) is actually analytic inside the fundamental 
domain. Therefore, we can conclude that the function f{u) must be zero function, i.e. 
f{u) = 0. This means that (I5.1UI) holds for I = N. 

Hence we have completed the proof of (15. 6p . 

Finally from the expression (15. 3p . we obtain the determinant representation of the par- 
tition function ZN{{ua}; {^i}] X; C) (I2.13P of the eight-vertex model with a non-diagonal 
reflection end and the DW boundary condition 

7 (I. \ ic\ \ /■\ TJ ^(-^12 + 2kr])a{Xi2 - 2kT] + 7]) ajuj + Ci) 

^NiiUahi^iii^i^ - II — 7T — : 1 \ /\ — ] — \~ 1 1 1 1 ~? — rr~^ — ^ 

IJl a{Xu + kr])a{Xi2 - kri + 7]) f^\f^\(^{Ui + ^i + V) 
11^=1 UiLi ^(^« ~ ^i)^iua + ^i + v) deW{{ua}; {^j}) 

Ua>f3 (T{Ua-Up)a{Ua+Up+v) Uk<l f^(^fc-6)^(^fe+6) ' 

(5.18) 

where the N x N matrix J\f{{ua}; {6}) is given by (15. 7p . 

6 Conclusions 

We have studied the partition function ZN{{ua}', X; () of the eight- vertex model with a 
non-diagonal reflection end defined by the K-matrix (12. 9p and the DW boundary condition. 
The DW boundary condition is specified by four boundary states (I2.38p - (l2.4ip . With the 
help of the F-basis provided by the Drinfeld twist for the eight-vertex SOS model, we obtain 
the complete symmetric expressions (15. 3p - (15.40 of the partition function. Such an explicit 
expression allows us to derive the recursive relation (15. 5p . Solving the recursive relation, 
we find that the partition function ZN{{ua}', {Ci}] X; () can be represented as a single deter- 
minant (I5.18p . contrary to the eight- vertex model without reflection end [131 [El [IS]. Such 
a single determinant representation will be crucial for constructing the explicit expressions 
of the Bethe vectors and further investigation of scalar productions of the open XYZ chain 
with non-diagonal boundary terms [52]. 
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